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Theoretical calculations, based on hybrid exchange density lunctional theory, are used to show that 
in graphene a periodic array of defects generates a ferromagnetic ground state at room temperature 
for unexpectedly large defect separations. This is demonstrated for defects that consist of a carbon 
vacancy in which two of the dangling bonds are saturated with H atoms. The magnetic coupling 
mechanism is analysed and found to be due to an instability in the n electron system with respect 
to a long-range spin polarisation characterised by alternation in the spin direction between adjacent 
carbon atoms. The disruption of the 7r-bonding opens a semiconducting gap at the Fermi edge. 
The size of the energy gap and the magnetic coupling strength are strong functions of the defect 
separation and can thus be controlled by varying the defect concentration. The position of the 
semiconducting energy gap and the electron effective mass are strongly spin-dependent and this 
is expected to result in a spin asymmetry in the transport properties of the system. A defective 
graphene sheet is therefore a very promising material with an in-built mechanism for tailoring the 
properties of the spintronic devices of the future. 

PACS numbers: 73.22.-f, 75.75.+a, 75.50.-y, 72.25.Dc, 73.90.+f, 75.30.Et, 75.30.Gw, 75.40.Mg, 75.50.Dd, 
75.50.Gg, 75.50.Kj, 75.50.Pp, 75.70.Ak 



I. INTRODUCTION 

The miniaturisation of silicon devices has delivered ex- 
ponentially growing performance for the past fifty years 
but quantum limits are now being reached as transistor 
gate lengths approach 5nm. Spintronics, which exploits 
the electron spin rather than the charge, will operate on 
a sub lnm length scale and thus facilitate another gen- 
eration of devices. 

The achievement of many of the promises of spintron- 
ics depends on the ready availability of a material that 
is a semiconductor, ferromagnetic at room temperature, 
and with a highly tunable band gap and magnetic cou- 
pling. The controlled synthesis of a robust room temper- 
ature, ferromagnetic semiconductor on a nanoscale has 
yet to be achieved. Recent reports of high temperature 
ferromagnetism in metal-free carbon materials have gen- 
erated a great deal of interest 3 . Proton-bombardment 
of graphite^ gives rise to a ferromagnetic phase which 
x-ray circular dichroism measurements^ have shown to 
originate from carbon sp electrons. Pyrolytic treatment 
of organic compounds^ also gives rise to a ferromagnetic 
matetieal. In both cases, the magnetic phase appears to 
be a minority phase in a non-magnetic host. The de- 
tailed composition and structure of the minority phase 
have not yet been determined but the concentrations of 
magnetic impurities have been measured sufficiently ac- 
curately to suggest that the magnetism is not the result 
of contamination^. 

The occurrence of high temperature ferromagnetism 
in purely sp-bonded materials is a major challenge to 



current theoretical understanding of magnetic interaction 
mechanisms. A chemical or structural defect in a carbon- 
based material leads to the creation of local magnetic mo- 
ments due to the presence of under- or over-coordinated 
atoms. The presence of local moments needs to be ac- 
companied by strong long-range coupling between them 
for the ferromagnetic order to survive thermal fluctua- 
tions at room temperature. In extended graphene sheets, 
local moment formation at point defects is now well es- 
tablished 8 . A significant step forward in the study of 
the magnetic coupling between local defects has been 
the demonstration that the presence of spin moments 
at the edges of graphene ribbons can lead to an insta- 
bility of the 7r-electron system, providing a mechanism 
for long range antiferromagnetic coupling which is robust 
with respect to competing instabilities involving charge 
ordering and geometric distortions ^ 10 ' 11 . More recently, 
Yazyev and Helmi2 have studied defective graphene us- 
ing density functional theory. They find a ferromagnetic 
ground state which they assume to be of the itinerant 
type due to the metallicity of their system. Vozmedi- 
ano et alJ£ studied the interaction between localised mo- 
ments originating from lattice defects at large distances 
and concluded that the transition temperature to the fer- 
romagnetic RKKY state is far below room temperature. 
Dugaev et alJS- predicted that if the spin-orbit interac- 
tion is taken into account, an energy gap opens at the 
Fermi level and the RKKY model is not applicable. A 
possibile source of high temperature ferromagnetism in 
carbon systems has been ascribed by Edwards et al^ 
to a limited screening of the interactions between itin- 
erant electrons in impurity sp-bands. The expected sup- 
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pression of the Stoner critical temperature (typical in the 
case of transition metal impurity d-bands) is considerably 
reduced due to a less effective screening of the electron 
interaction. In the present case we show a completely 
different scenario, i.e. that high temperature ferromag- 
netism in a carbon system can be entirely supported by 
an insulating state of the electronic structure. 

The current work demonstrates that local spin mo- 
ments in defective graphene interact strongly over long 
distances and may therefore lead to a ferromagnetic semi- 
conducting state at room temperature. The strong, long- 
range magnetic coupling necessary for magnetic order- 
ing at room temperature in this semiconducting material 
is provided via a spin alternation mechanism that de- 
pends on the presence of partially filled tt orbitals in the 
sp2-bonded, bipartite lattice of graphene. Moreover, it 
is shown that the spin-dependent band gap and mag- 
netic coupling strength can be varied through control of 
the defect concentration, opening up the possibility of 
highly tunable graphene-based spintronics nanodcviccs. 
This study, therefore, reinforces the case, which has been 
growing stronger and stronger since the ground-breaking 
practical realisation of isolated graphene sheets 2 -, for con- 
sidering graphene amongst the most promising materials 
for building the electronic and spintronic devices of the 
future. 



parallel chains running along the horizontal direction and 
the chains are 25 A apart. The second structure is an 
hexagonal cell used to quantify inter-chain effects. The 
third structure is a triangular superlattice used for rep- 
resenting the band structure in a form compatible with 
that of graphene (Fig. 4) . All structures considered were 
initially relaxed using a force-field model (GULP—) in 
which the carbon-carbon and carbon-hydrogen interac- 
tions are described via the Brenner bond order poten- 
tial^. In the case of the calculation of the stabilisation 
energy AE (see below and Fig. 2), full quantum me- 
chanical relaxation was carried out for the structure with 
L=7.5A (see below) and resulted in a small increase in 
AE (see Fig. 2) of 25 meV. As the additional relaxation 
is confined to the region of the defect, the energy shift 
is expected to be depend weakly on the defect separa- 
tion especially at large distances and thus full quantum 
mechanical relaxation was not required for larger struc- 
tures. The magnetic stablisation energy reported in Fig. 
2, and the transition temperature based on it, is there- 
fore slightly underestimated. In the analysis of the elec- 
tronic band structure shown in Fig. 4, where defects 
are displayed in a triangular lattice and are 20 A apart, 
we performed a full quantum-mechanical optimisation in 
order to quantify the magnitude and decay law of the 
semiconducting gap. 



II. COMPUTATIONAL DETAILS 

Spin localisation at defects and long range spin 
coupling are studied here using all electron, two- 
dimensionally periodic, hybrid exchange density func- 
tional theory (B3LY P 17 i 18 ' 19 ) which significantly extends 
the reliability of the widely used local and gradient cor- 
rected approximations to density functional theory in 
strongly interacting system a 20 ' 21 i 22 . In the CRYSTAL 
package, used for this study, the crystalline wavefunc- 
tions are expanded as a linear combination of atom cen- 
tred Gaussian orbitals (LCAO) with s, p, d, or / sym- 
metry. The calculations reported here are all-electron, 
i.e., with no shape approximation to the ionic potential 
or electron charge density. Basis sets of double valence 
quality (6-2IG* for C and 6-31G* for H) are used. A 
reciprocal space sampling on a Monkhorst-Pack grid of 
shrinking factor equal to 6 is adopted after finding it to 
be sufficient to converge the total energy to within 1CU 4 
eV per unit cell. 

In the present work we consider three types of struc- 
ture, characterised by different arrangements of defects. 
The first structure contains defects arranged in pairs 
within a rectangular unit cell that is 25 A wide and of 
length 2L, which is convenient for studying inter-defect 
coupling. The system under study is therefore infinitely 
periodic in the 2 dimensions parallel to the graphene 
plane and non-periodic in the direction perpendicular to 
the plane. In Fig. 3 the unit cell containing a pair of de- 
fects is shown. The defects are therefore arranged along 



III. RESULTS AND DISCUSSION 

A number of point defects are possible. Here we choose 
to study a simple carbon vacancy in which two of the 
three dangling carbon bonds are saturated with H-atoms. 
This defect has been proposed previously, studied exten- 
sively*, and is displayed in Fig. 1. Within a valence 




FIG. 1: The structure of the defect considered in the present 
study. 

bond picture there are two unpaired electrons in the de- 
fect as the breaking of three 7r-bonds generates 1/3 of an 
unpaired electron on each of the C-atoms neighbouring 
the vacancy and the breaking of the cr-bond also gener- 
ates an unpaired electron. Hund's rule coupling within 
the uncoordinated C-atom and the spin-alternation rule 
operating between the 7r-electrons ensures that these 
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FIG. 2: The dependance of the magnetic stabilization energy 
AE on the distance between defects L. 



spins are aligned and thus the overall magnetic moment 
on the defect is 2/xb (Ms is the Bohr magneton). — In the 
current study this defect is used as an example of the 
many possible structures at which electron localisation 
might occur in a graphene sheet. The conclusions drawn 
below depend on the arrangement of defects and on their 
ability to generate a local magnetic moment through an 
instability of the electronic system, but are not expected 
to depend on the structural and chemical details of the 
defects. 

The interaction of the localised moments is quanti- 
fied by computing the quantity AE (stabilization energy) 
which is defined as the difference between the total elec- 
tronic energy of the antialigned and aligned arrangements 
of the defects' spins. The decay of the stablisation en- 
ergy with L is approximately AE = 4208L -143 meV 
and it is shown in Fig. 2. The calculated coupling en- 
ergy can be related to the Curie temperature (T c ), at 
which the ferromagnetic order is lost due to thermal fluc- 
tuations, through an effective Heisenberg model with an 
interaction J = AE/2 per pair of spins. According to the 
Mermin- Wagner theorem 24 , a two dimensional isotropic 
Heisenberg ferromagnet is unstable with respect to ther- 
mal fluctuations at any finite temperature and therefore 
does not support long range order. However such effects 
are seen to reduce significantly when a small coupling 
along the third direction is introduced. Recent Quantum 
Monte Carlo (QMC) studies on quasi two dimensional 
Heisenberg ferromagnet a 25 ' 26 have shown that the intro- 
duction of an interlayer coupling 0.1% of that in-plane 
is enough to restore a finite value of T c . In particular, 
for a square lattice of spin-1 T c = JS(S + l)QAA/K B 
(see Fig. 3 in Ref. 26 ) which, compared to a mean-field 
expression T c = JS(S + 1)z/3Kb (z is the coordinaton 
number), differs by a factor of about 3. On the right-hand 
side of the graph in Fig. 2, the transition temperature 



T c for a square lattice arrangement of defects within the 
QMC formula above is shown for different interdefects 
distances. For this arrangement of the defects, the ferro- 
magnetic state is lower in energy over the whole range of 
distances considered and is predicted to be stable at room 
temperature up to very large defect separations (Fig. 2). 

The present system should be considered as a model 
system for the study of more complex and realistic 
cases such as defective graphite^ or multilayer defective 
graphene or defective graphene on a substrate where a 
small but non zero exchange interaction in the third di- 
rection is usually found. The ferromagnetic long range 
order at room temperature is in agreement with recent 
experiments on bombarded graphite where a room tem- 
perature ferromagnetic ground state is found and the av- 
erage inter-defect distance was determined to be about 1 
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Effects due to interchain interaction expected at dis- 
tances L close to the interchain fixed distance of 25 A are 
seen not to change the picture qualitatively. To confirm 
this, we independently considered a series of calculation 
of AE for a strictly 2-dimensional arrangement of defects 
(a hexagonal superlattice of defects) for which the quali- 
tative nature of the long ranged interaction is unchanged 
and the power of the decay law is found to be 1.7. 

The bipartite nature of the lattice and the presence of 
the partially filled it system of electrons in the sp 2 net- 
work of atoms are essential for establishing the strong, 
long-range coupling mechanism that gives magnetic or- 
dering at high temperature. This mechanism is the result 
of local exchange repulsion, which ensures that the par- 
tially filled 7r orbital of each carbon atom is coupled in 
an anti-parallel fashion to that on its nearest neighbour. 
The operation of this mechanism to produce a long range 
spin polarisation within the sheet is evident in the spin 
density displayed in Fig. 3 and has been observed previ- 
ously in 7r bonded carbon networks 28 . It is notable that 

this mechanism extends the validity of Lieb's theorem be- 
yond a nearest neighbour model 29 . The predicted ferro- 
magnetic nature of the overall interaction is also entirely 
consistent with previous studies of interactions between 
extended defects in graphene ribbons2J£. It is clear from 
the spin alternation mechanism that the ferromagnetic 
coupling between the localised spin moments depends 
on the defects being arranged on the same sublattice of 
the graphene sheet. A regular arrangement on opposite 
sublattices would lead to antifcrromagnetic coupling, as 
seen by Yazyev et al. at short separation^ 2 ., and a ran- 
dom arrangement of defects would be expected to lead 
to a ground state of smaller spontaneous magnetisation, 
while retaining the strong local coupling. This predic- 
tion agrees with the current experimental observation by 
Esquinazi et. al££Z- (i.e. high transition temperature 
and small total magnetisation), where the defect density 
is controlled, to an extent, via the defect implantation 
technique but the defects are most probably distributed 
randomly. The production of a graphene sheet and con- 
trol of its defect density is not unrealistic using current 



4 




FIG. 3: (color on line) Isovalue surfaces of the spin density of 
the graphene sheet with a defect separation of 20 A. The red 
(at 0.021 ps/A) and blue (at -0.021 ps/A)) isosurfaces rep- 
resent the majority and minority spin densities, respectively. 
The majority spin density within the defect is concentrated on 
the unsaturated carbon atoms and their second neighbours. 
The operation of the spin alternation rule in the spin polarised 
lattice is clearly visible. 



techniques 2 -^ but methods for controlling the defect ar- 
rangement have not yet been established. 

To analyse the semiconducting properties of the elec- 
tronic structure of the system considered in the present 
study, the electronic band structure is show in Fig. 4 for 
the ferromagnetic state of a triangular array of defects 
at a separation of 20 A . This choice of array facilitates 
the comparison with that electronic structure of perfect 
graphene. The previously analysed structures (parallel 
chains and hexagonal arrangement of defects) show the 
same qualitative features as for the triangular array. 

The presence of the defects and the magnetic ordering 
break the symmetry of the graphene 7r system, opening 
band gaps of 0.51 and 0.55 eV in the majority and mi- 
nority spin band structures, respectively. The energy gap 
of the majority states is shifted upwards by about 0.20 
eV with respect to the energy gap of the minority states, 
and the shape of the bands around the Fermi energy is 
markedly spin dependent. A strong spin asymmetry in 
the conductivity of the sheet is therefore to be expected. 
The size of the band gap is reduced to 0.1 eV when adopt- 
ing the PBE generalised gradient approximation^. 

Calculations at a variety of defect separations estab- 
lish that the band gap scales as L~ 2 . This decay law is 
consistent with that recently calculated for graphene rib- 
bons where the edges are line defects at which the spin 
moments localise and the band gap scales as L _1 with 
the ribbon width. 

Recently Yazyev et al^ reported GGA calculations 
that produced ferromagnetism in graphene induced by 
the presence of defects but the magnetic coupling mecha- 
nism proposed differs strongly from that described in the 
current work. In Rcf. 12 a metallic state is obtained and 
the ferromagnetic order is therefore assumed to be sup- 
ported by itinerant electrons. In the present study the 
ground state is semiconducting and the magnetic cou- 




FIG. 4: (color on line) The electronic energy bands of the 
majority (a) and minority (b) spin states of the ferromag- 
netic ground state in a defective graphene sheet with a defect 
separation of 20 A plotted with respect to the Fermi energy 
(E_f). The rather flat impurity bands near the Fermi energy, 
which are associated with the states localised at the defect, 
are indicated in red. 



pling is clearly not due to itinerant electrons. For this 
reason we adopt an effective Heisenberg model in order 
to assess the temperature related properties of the sys- 
tem rather than the Stoner model of Reffi 2 ". The ex- 
planation for this different behaviour does not reside in 
the morphology or in the concentration of defects nor 
in the treatment of electronic exchange and correlation 
(the band gap is still present when using the gradient cor- 
rected functional PBE as adopted in RefJ^) but is due to 
the arrangement of defects. In Ref^, the authors study 
the same H-chemisorption defect as in Refj^ 2 -, but con- 
sider a different arrangement of defects (specifically, the 
unit cell in Ref^ has zig-zag borders whereas armchair 
borders are used in Ref^). This results in a semicon- 
ducting ground state as found here. In Rcf. 30 the mag- 
netic coupling between defects was not analysed and this 
is the aim of the present study. 

The opening of a gap obtained in this study is an effect 
of much larger magnitude than the opening due to spin- 
orbit coupling described in RefJ^. The effect of the spin- 
orbit coupling is not included in our study since it is 
expected to be negligible in carbon systems 3 ^,. 

IV. CONCLUSIONS 

In conclusion the ground state of a two dimensional 
graphene sheet containing a periodic array of point de- 
fects has been shown to be both ferromagnetic at room 
temperature and semiconducting for defect separations 
up to 20 A and thus for defect concentrations as low 
as 10 13 /cm 2 . The energy gap and magnetic coupling de- 
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pend strongly on defect concentration. We conclude that 
a doped or defective graphcne sheet is a very promising 
material with an in built mechanism for tailoring prop- 
erties for a variety of spintronics applications. 
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